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Abstract: Using a “risk network” to describe the contact patterns (including needle sharing and 
sex behavior) of buddy relationships in injecting drug users (IDUs), a stochastic model 
is studied about the spreading of HIV-1 in IDUs. The basic reproduction number and 
the final size of the epidemic are obtained. In particular, the model has been applied to 
Xichang city. Through simulations, the effect of the injecting network size (and also the 
sex network size) in the spreading of HIV-1 in IDUs are studied. The model shows that 
10% increasing of the final size of the epidemic need only one triple of needle sharing added 
for one IDU. The results give us very important hint: reduce the network sizes, especially 
the needle sharing network size should be a very important intervention measure to control 
the spread of HIV-1 in IDUs. 
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1 = Introduction 


HIV-1 is an infectious pathogen that is now largely spread through close personal contract. 
A “risk network” is comprised of those people among whom HIV-1 risk behaviors occur. Risk 
networks thus act as conduits of infection. The role of risk networks as conduits of HIV-1 infec- 
tion is illustrated by studies!'-3! that have focused on IDUs in Sichuan province and Xinjiang 
province, China. Qualitative data were collected and were used to construct sociograms of the 
risk networks of drug injectors. Authors found that the size of risk network of sharing injection 
equipment among IDUs has positive association with HIV-1 infection. We use the term “risk 
network” to describe the contact patterns (including needle sharing and sex behavior) of buddy 
relationships that are stable over a longer time. So the problem here is how to find a graph 
that captures the social formations in drug users in a satisfactory way. 

IDUs in China normally have two risk social relationships which are suited for the spread 
of HIV-1: one is the “buddy relationship” (needle sharing network), and each group size is not 
very big (normally include about 2 to 3 people). The other is the “sexual relationships”. We 
will use an undirected graph to describe the spread of the epidemic. Given a set of N labeled 
vertices. The risk network will be constructed in the following way: 
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1) We add edges in all vertices (v;,v;) with probability rq, and the subgraph that is 
generated is denoted by G4; 

2) We add edges in all possible triples of vertices (vi, vj, Vk) hjk =1, N, tF GFK 
(call them triangles) with probability fa, and the subgraph that is generated is denoted by G4, 
Both G? and G® represent the needle sharing network; 

3) We add edges randomly in all vertices (v;, vj), i,j = 1, , N, 044 with probability rs, 
and the subgraph that is generated is denoted by G*. Let it represent the “sexual relationship” 
structure, 

So this construction procedure generates a graph G that can be written as the union of three 
independent subgraphs G4, G4 and G5. Also, we call rs, fa and rg as network probabilities. 

Given a graph G generated as above. The average coordination number z for a given vertex 
v in G is 

= » ri(N —1)+ NP. 
i=s.d 


Since no people can have infinite friends, we put 
re=M/(N-1), i=sd, fa=da/(CY”) 


for some A; > 0 (i = s, d) and da > 0 to scale the edge probability in G. Here Aa, As and dg are 
the average numbers of two-people fellowship of needle sharing. sexual fellowship and triangle 
fellowship of needle sharing of vertex v respectively. 

In the following we discuss what will happen if we let HIV spread along graph G. We 
consider a closed population, consisting of N individuals. Let n denote the initial number 
of susceptible individuals and m (which is very small compare with n) be the initial number 
of infectious individuals. All individuals in the network are labeled as either susceptible (5), 


infected (I) or removed (R). 


2 The basic reproduction number and the final size of the epidemic 


First we assume that addicts sequentially use the same drug-injection equipment in G4, 
Any uninfected addict who uses infectious injection equipment is considered to be exposed 
to HIV-1 of course a given HIV-1 infective can only infect his buddies (e.g., his sex partners 
by probability ps or his drug injecting parters by probability a and pq). Also the infection 
probabilities should be thought of as a product of the contact probability and the probability 
of a disease transmission in case of a contact, thus giving the probability of a contact resulting 
in a new infective. 

Now provided that the population is large, contacted individuals in the epidemic process 
are susceptible with high probability in the beginning of the time course. Thus the initial 
stage in the generation process of infective is well approximated by a branching process with 
reproduction mean equal to the average number of new cases generated by a given infective. We 
denote it as the basic reproduction number Ro. This is a critical parameter indicating whether 
or not a large outbreak is possible. If Ro < 1 then the asymptotic probability of a large outbreak 
is zero, while if Ro > 1, then there is, asymptotically, a strictly positive probability of a large 
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outbreak. By a large outbreak we mean an outbreak such that the size of the outbreak is of 
the same order as the whole population. 

For G = G°UG4UG". Let Rs, Rg and Ry denote the expected numbers of secondary cases 
generated by an infected individual 7 in G°, G4 and G4 respectively. Then Ro = Ry + Ra + Ra. 

Let X° and X? denote the degree of vertex i in G° and G4, respectively, and X‘ is the 
number of triangles in G4 involving vertex i. 

Since each of these individuals is infected by i with probability ps in G, for sex active people, 
hence Rs = psE[X*] = pg. Similar can get that Ra = paE[X4 = pada. 

About Ra, it should be the product of two expectations: one expectation is of the number 
of triangles in G4 involving i, that is, E|X gos da; the other expectation is of the number 
of secondary cases generated by į in a fixed triangle, we denote it as E[Ý¢2]. All the possible 
values of Y¢ should be 0, 1 or 2. It is easy to get that the probability of only one new infective 
is 2pa(1 — pa)? since a is the transfer probability by the 3-clique sharing injecting equipment 
group. In case of two new infectives, there are two possibilities: the two new infectives are 
infected directly by i, which probability is Di; or the infection takes place in two steps, which 
probability is 2(1 — pa)p2. So 


E(Y4] = 2pa(1 — pa)? + 2[54 +2(1 — pa)54], Ra = E[X¢E[Y4] = 2\a( Da + 4 — B3). 


Let R? be the total reproduction number of sharing injecting equipment groups including 
these 2-cliques and these 3-cliques, that is R? = Ra + R4, then we have 


Ro = Rs + RY = Agps + 2alBa + P2 ~ PS) + Aapa. 


In the following, we study the final size of the epidemic. Suppose the initial infected individ- 
uals number No is very small. Define Ts, 7g and Tq as the proportion of the initially susceptible 
individuals in G*, G4 and G4, respectively, that is ultimately infected. Define 7 as the propor- 
tion of the whole initially susceptible individuals that is ultimately infected. Asymptotically the 
probability that a given individual is infected should equal to the proportion of the population 
that is ultimately infected in the epidemic. 

Now consider a given susceptible individual i. Let A denote the event that i is infected. Let 
B, C and D denote the event that i is infected through G*, G4, or G4, respectively. Then 





P(A) = P(B)P(C)P(D), Ts = P(B), a= P(C), ta = P(D), 
T=T3+Ta+Ta — TsTa — TsTA — TaTa + TsTaTa- (1) 
Let B; denote the event that 7 is infected from a given neighbor j in G°. Then 
P(B;) = Taps + (1—75)0= Tops, P(B) = E[P(B;)*"] = E[(1 — reps)*’]. 


Take the right-hand side in the above equality as the probability generating function of X* 
evaluated at 1—7,p,, then Ts = 1—e77sPsàs, Using similar way we can get that rg = 1—e7~T4PaAa, 


Now we consider P(C). Since we suppose that the number of triangles in G4 involving 


vertex i is Xt, so 


PO) = E| P(N ANTk)] = ALPE], v=  Q) 
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where Cy means that 7 is escapes infection in a fix triangle V = {i, j, k}. 
Let A denote the event that j and k will be both infected as soon as one of them is infected. 
Then 


P(Cv) = P(Cy | A) P(A) + P(Cy | A°) P(A‘). (3) 


Now we have three possible situations: both k and j are not infected; one of j and k is 
infected; and both k and j are infected. According to our definition of Ta, the asymptotic 
probability for an individual to escape infection is equal to 1 — 7y. So from (3) we can get 


P(Cy) = pall — 7a)? + 27a(1 — 74) ~ 4)? + 73.0 — pa)?) 
+(1 = pa) [0 — Fa)? + 27a(1 ~ F4)(1 — Ba) + — pa)? 


1 ~ [2pata(2 — Fa) + ByFa(2 — 374) — 23Fa(1 — Fa). 


II 


Looking at the right-hand side in equality (2) as the probability generating function of X4 
evaluated at P(Cy), also X? is asymptotically Poisson (a) distributed, then 


CO) = edal2Bata(2—Fa) +53 Fa (2—3F4) ~ 23 Fa(1—-Fa)] 
P(C) =e ; 


So 


Fy = 1 — e`% l2Pa?a (2-a) +pa3a(2—374) -247a (1~Fa)] 
and equation (1) holds too, that is 


T=TtTa+ Ta — TsTa — TsTd — TiTa + TsTdTd- 


3 Applying to Xichang city 
To get the infection probabilities ps, pa and pa, we define 
Ps = hsns(1 -— ests), Pa = haña(l — eaca), pa = hana(l — eaca), 


in which, hs is the transmissibility by sex connection, and ns is the sexual frequency between 
the infected individual and one of his sex partner. The transmission of HIV-1 also depends 
on condom efficacy, €s, and the frequency of condom use, cs, between the infected individual 
and his sex partners. We use hg represent the transmissibility that an addict becomes infected 
given exposure to HIV-1. It is the infectivity of HIV-1 via shared injection equipment. fig and 
Na denote the sharing frequency of injection equipment between the infected individual and one 
of his friends in G4 and G4 respectively. Also we assume that an addict cleans the needle with 
probability cą prior to use, and parameter eg means the efficacy of the cleaning. According to 
literatures!*-9), we choose the following parameters values as that in F igure ł and Figure 2 for 
a given vertex v in G. 

First we discuss how the injecting network size (and also the sex network size) effect the 
outcomes of HIV spreading. Let Às, Ag and Xa denote the sizes of the sex network and the 
injecting network. The relationship between the reproduction number Ro and the sex network 
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size À, can be seen in Figure 1. This figure shows clearly that the reproduction numbers Rs and 
Ro increase quite fast along with the increasing of the sex network’s dimension. The star points 
are the reproduction numbers of the parameters in Figure 1 and A, = 3. Similar trend curves 
can be got for the relationship between the reproduction number Ro and the needle sharing 
network size of IDUs, Ag. 




















Figure 1: R, and Ro increase fast along with the increasing of the sex network size As. Other 
parameters are chosen as da = 2, Ad = 3, hs = 0.05, ha = 0.3, ns = 8, es = 0.7, cs = 0.09, na = 1.3395, 
ña = 0.893, ca = 0.04, eg = 0.6. 
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Figure 2: Tq and T increase very fast along with the increasing of the risk network size da (i.e., 
xa"): Other parameters are chosen as As = 3, àa = 3, hs = 0.05, ha = 0.3, ns = 8, es = 0.7, cs = 0.09, 
Na = 1.3395, ña = 0.893, ca = 0.04, ea = 0.6. 
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The importance of the size of the risk network in the spreading of HIV-1 in IDUs, especially 
for the number of sex partners and the needle sharing triple, can be seen from simulations. 
The size of needle sharing network (which includes these triples) contributes much more to 
the spreading of HIV-1 in IDUs than that of the sexual network. Figure 2 shows that, 10% 
increasing of the final size of the epidemic need only one triple of needle sharing added for one 


IDU. 
So we conclude that, reduce both of the network sizes, especially the needle sharing network 
size should be a very important method to control the spread of HIV-1 in IDUs in China. 
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BF fe he ee ba A BE HIV BPE 
ih, PR X 


(LAFF AR, ie 200444) 


A E WAU RBAR PIERR (AHR AMERA) H “SCM KRMR” , RIMAT -AR 
F HIV ELAR EOL, BET RAE AR RS el. EO GS a 
E, RAVI T HARA PERN PER ARS ER. UL, ENE E 
WIM—t “= A/D” WR, BA SB 10% KBR A. ARER: RZD FE RP UA, 
TUE ASAI a, HFN HIV eR EAD 

Hi: HIV/AIDS; ERAR: EREM IDUs 





